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We study the quantum phase transition of the Dicke model in the classical oscillator limit, where 
it occurs already for finite spin length. In contrast to the classical spin limit, for which spin-oscillator 
entanglement diverges at the transition, entanglement in the classical oscillator limit remains small. 
We derive the quantum phase transition with identical critical behavior in the two classical limits 
and explain the differences with respect to quantum fluctuations around the mean-field ground state 
through an effective model for the oscillator degrees of freedom. With numerical data for the full 
quantum model we study convergence to the classical limits. We contrast the classical oscillator limit 
with the dual limit of a high frequency oscillator, where the spin degrees of freedom are described 
by the Lipkin-Meshkov-Glick model. An alternative limit can be defined for the Rabi case of spin 
length one-half, in which spin frequency renormalization replaces the quantum phase transition. 
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I. INTRODUCTION 

For a system of a single spin coupled to a quantum har- 
monic oscillator a quantum phase transition (QPT) [1 
can take place only in the classical limit of one of the 
two components, i.e. in the limit of infinite spin length 
or zero oscillator frequency. Prior to the respective clas- 
sical limit, the spin-oscillator system admits no phase 
transition since symmetry breaking states can always be 
combined in a linear superposition that restores the sym- 
metry and reduces the energy further. In the classical 
limit phase transitions become possible because different 
classical states have zero overlap, which circumvents the 
previous argument against symmetry breaking. 

This type of QPT is realized in the Dicke model [2] 

H = AJ Z + Qala + A(a f + a)J x . (1) 

It describes an ensemble of 2j two-level atoms with tran- 
sition frequency A as a pseudo-spin of length j (using 
spin operators J x / Z ). The atoms are coupled to a single 
cavity mode of the photon field with frequency (using 
bosonic operators a^'). The Hamiltonian in Eq. ([I]) is in- 
variant under the replacement J x > — J x , a t— > —a. This 
symmetry is broken in a phase transition, and the spin 
expectation value ( J x ) serves as the order parameter. 

In the classical spin (CS) limit j — > oo the Dicke model 
features a thermodynamic phase transition from a high 
temperature state with ( J x ) = (a) = to a superradiant 
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state with a finite cavity field ((a) ^ 0) and macroscopic 
atomic excitation ((J x ) ^ 0) at low temperatures [3H5]. 
The thermodynamic phase transition is complemented 
at zero temperature by a QPT from the zero field to 
the superradiant state at a critical atom-field (i.e. spin- 
oscillator) coupling A c . 

The driving mechanism behind the QPT is the criti- 
cal behavior of a classical energy functional for the spin, 
which is obtained after integrating out the quantum- 
mechanical oscillator. Strictly in the j = oo-limit the 
ground state is a mean-field (MF) product state of a spin 
and oscillator coherent state. The order parameter ( J x ) 
and the corresponding susceptibility x, which character- 
ize the critical behavior, converge to the classical results 
if the CS limit is approached from j < oo. 

Modifications of the classical picture arise from quan- 
tum corrections of order 1/j to the MF ground state [SJ 
7 . Spin and oscillator variances diverge and signal the 
breakdown of the classical limit in the vicinity of the 
QPT. A characteristic feature is the criticality of spin- 
oscillator entanglement [8HTU]. which is related to the 
vanishing excitation gap at the QPT and found for many 
different models in the CS limit [TTj . 

In this paper we address a QPT in the different clas- 
sical limit £1 — > 0, the classical oscillator (CO) limit. In 
contrast to the CS limit, a QPT transition occurs here 
already at finite spin length j. The critical behavior is 
identical to the CS limit, since both limits realize the 
same MF transition. Quantum corrections to the MF 
ground state are different. In the CO limit spin fluctua- 
tions are suppressed because of the large spin frequency. 
Therefore, the spin variance and the spin-oscillator en- 
tanglement remain small in the vicinity of the QPT. The 
entanglement entropy is bounded by In 2 independently 
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FIG. 1. (Color online) Diagram of the qualitative behavior 
of the Dicke model in dependence on the spin length j and 
the spin/oscillator frequency ratio fl/A. The bold axes cor- 
respond to the QPT in the CS limit 1/j = and the CO 
limit n/A = 0. The top edge of the square corresponds to 
the FO limit fi/A — > oo, where the LMG model describes 
the spin. On the right side of the square we find the Rabi 
model (j = 1/2), with renormalization of the effective spin 
frequency in the limit fi/A — > oo (upper right corner). 



of j. The CS and CO limits thus give rise to QPTs with 
identical critical behavior that are distinguished through 
the criticality versus non-criticality of entanglement. 

An overview of the different limits in the Dicke model 
is given in Fig. [T] The paper is organized according to 
this diagram. We first derive in Sec. [TT| the QPT in the 
CS and the CO limit from MF theory, which becomes 
exact in the two limits. Quantum corrections in the CO 
limit are discussed in Sec. Mil with an effective bosonic 
model for the oscillator degree of freedom, and contrasted 
with the behavior in the CS limit known from litera- 
ture. The (non-) criticality of entanglement is addressed 
in Sec. |IV| We complement the CO limit with the fast 
oscillator (FO) limit fl/A — >• oo in Sec.|v] Similar to the 
considerations for the CO limit, the large oscillator fre- 
quency leads to the suppression of oscillator fluctuations. 
One obtains the Lipkin-Meshkov-Glick (LMG) model for 
the spin degree of freedom, but no QPT occurs unless 
we let again j — > oo. A related FO limit that is pecu- 
liar to the Rabi case j = 1/2 leads to renormalization 
of the spin frequency instead of a QPT. The appendices 
summarize the solution of the effective bosonic models 
for the CO and FO limit (App. [A| ), and the definition 
of coherent states (App. [B]) and of the rotation invariant 
spin variance (App. O). 



II. CLASSICAL LIMITS AND THE QUANTUM 
PHASE TRANSITION 



As noted in the introduction, the Dicke Hamiltonian 
from Eq. ([I]) is invariant under the symmetry transfor- 
mation 



n = e 



N E = a)a + J z +j , 



(2) 



which corresponds to the simultaneous replacement of 
a i y —a and J x n- — J x . We have n _1 .HTI = H or 
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FIG. 2. (Color online) Mean-field QPT in the Dicke model. 
Left panel: Energy functional E(9) from Eq. Q below (k < 
1), above (k > 1) and at (k = 1, dashed curve) the QPT. 
Right panel: Order parameter ( J x ) (red) and susceptibility % 
(black) of the mean-field QPT given in Eqs. flu], |TT|. 



[H, II] = 0, such that the eigenstates of H can be clas- 
sified by the eigenvalues ±1 of the parity operator II. 
For positive A, the ground state of the Dicke model has 
positive parity. 

The QPT breaks parity symmetry, with a finite order 
parameter (J x ) ^ above a critical coupling. To study 
the convergence to the classical limits it is useful to break 
the parity symmetry explicitly, using the Hamiltonian 

H £ =H~ e,J x (3) 

that includes a symmetry breaking field eJ x . 

A. Mean-field theory of the QPT 

Let us first discuss the QPT in the MF picture. The 
MF ansatz for the ground state wave function 



MF/ 



a) 



(4) 



is a product of a spin coherent state \0) and a boson 
coherent state \a) (cf. App. |Bj. This state has energy 

E(e,a) = {tPMFm M F) 

2 W 

= — j A cos 9 + fla + 2j\a sin 9 . 
Minimization with respect to a results in 



a = — — sin ( 



(6) 



which inserted into Eq. (|5| gives the energy functional 

E{9) = - j A (cos 9+^ sin 2 6>) . (7) 

We here introduce the dimensionless coupling constant 

2jA 2 



(8) 



We will assume k > 0, which corresponds to A > 0. 

The energy functional E(9) describes a second-order 
MF transition at k = 1 (cf. Fig. [2). The minima of E{9) 
are given by 
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if k < 1 
if k > 1 



(9) 
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which leads to the expression 
'0 



(Jx 



j sin t 



±3\ 1 



if k < 1 , 
if k > 1 



(10) 



for the order parameter (J x ). We can also calculate the 
susceptibility x using the Hamiltonian H e from Eq. ([3]), 
and find 



1 



X = jA lim 



d(J x ) 
6 de 



1 - K 
1 



k(k 2 — 1) 



if k < 1 , 



if k > 1 



(11) 



In contrast to the prediction of MF theory, the argu- 
ment given in the introduction shows that a QPT cannot 
exist in the fully quantum-mechanical Dicke model for fi- 
nite j, ft/ A. The QPT only becomes possible if the two 
degenerate MF states for positive/negative 9, a have zero 
overlap. This is can be achieved either if (0\—9) = in 
the CS limit, or if (a\ — a) = in the CO limit. 



B. QPT in the classical spin limit 
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FIG. 3. (Color online) QPT in the CS limit. Order param- 
eter (J x ) (left panel) and susceptibility \ (right panel) as a 
function of k, for various values of j and fixed fi/A = 1. The 
dashed curves give the MF result from Eqs. ( 10 1, (111. 
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FIG. 4. (Color online) QPT in the CO limit. Order param- 
eter (J x ) (left panel) and susceptibility % (right panel) as a 
function of k, for various values of fi/A and fixe d j — 5. The 
dashed curves give the MF result from Eqs. ( 10 1, (111. 



In the CS limit j — > oo spin coherent states form an 
orthonormal basis of the spin Hilbert space (see, e.g., 
Ref. [H]). Therefore, the ground state wave function has 
the form \ipcs) = \9) ® IV'bos ), w ith a real spin coherent 



state \9) as defined in Eq. (JB4J) in App. |Bj Note that 
for 9 7^ the overlap (9\— 9) = cos 2 - 5 9 goes to zero for 
j — » oo, which allows for the QPT. 

The bosonic part \ipbos) of the wave function, which 
has to be determined through minimization of the energy 
(tpcs\H\tpcs) , is the ground state of the effective bosonic 
Hamiltonian 



H osc (9) = f2a f a + Xj sin 9(a + a)) , 



(12) 



which is parameterized by the classical spin angle 9. 

H osc (9) is the Hamiltonian of an oscillator with a con- 
stant force oc Xj sin 9. The ground state of this Hamilto- 
nian is a boson coherent state |^bos) = with a given 
by Eq. ([6]). We thus recover the MF wave function from 
Eq. Q in the CS limit, hence also the entire QPT. 



C. QPT in the classical oscillator limit 



of the wave function is the ground state of the effective 
spin Hamiltonian 



H spin (a) = AJ Z + 2XaJ x 



(13) 



parameterized by the classical oscillator displacement a. 

H S pi Yl (a) is the Hamiltonian of a spin in a magnetic 
field B = (2Aa, 0, A), with a coherent spin state |V>spin) = 
\9) as the ground state. 9 and a are related through 
Eq. (6{. Again, we recover the MF wave function from 
Eq. (4J), and therefore also the entire QPT in the CO 
limit. 

Note that the argument for the CS and CO limit are 
dual to each other: In the CS limit we first observe that 
the spin has to be in a coherent ("classical") state and 
deduce the oscillator coherent state from the particular 
effective Hamiltonian H osc (9) for the quantum mechan- 
ical oscillator. In the CO limit, we start from an oscil- 
lator coherent state and obtain the spin coherent state 
again only because of the particular form of the effective 
Hamiltonian H sp - ln (a) for the quantum spin. 



According to Eqs. p|, ^ the parameter a in the MF 
ground state in Eq. Q scales as For fl — > 

the overlap (a>\— a) = exp(— 2a 2 ) goes to zero for k > 
1. Because of this a QPT in the CO limit is possible 
independently of the spin length j. 

Since the overlap of different real coherent states \a) 
is zero in the CO limit, the ground state wave function 
has the form \ipco) = l^spin) ® \a). The spin part |^ sp in) 



D. Convergence to the QPT 

Given the QPT in the CS and CO limit, we expect 
that the order parameter (J x ) and the susceptibility x 
converge to the MF values from Eqs. ([10]), ([IT]) if the 
classical limits are approached from the quantum regime 
j < oo, fl/A > 0, i.e. if one moves in Fig. [l] from the 
interior of the square towards one of the two bold axes. 



4 



To observe convergence we must use a small symmetry 
breaking field eJ x as in Eq. ^ to select one of the two 
possible cases (J x ) ^ in the broken symmetry phase 
K > 1. Without the additional field convergence cannot 
be observed because parity symmetry implies that ( J x ) — 
in the quantum regime. Apart f rom the calculation of 
the entanglement entropy in Sec. 
throughout the paper. 



IV 



we use e 



In Figs. [3j [3] we show (J x ) and % from a numerical cal- 
culation of the ground state of the Dicke model using the 
Lanczos technique [13 . Up to 10 3 bosons are kept in the 
calculations to ensure a negligible error from the trun- 
cation of the infinite-dimensional bosonic Hilbert space. 
The spin part is not truncated, and the numerical data 
are accurate on the level of machine precision. 

In the figures we start from the curve for j = 5, 51/ A = 
1, far away from the classical limits. In Figj3]we increase 
j to approach the CS limit, and in Fig. |4pwe decrease 
51/ A to approach the CO limit. In both situations we 
observe convergence of (J x ), x to the MF values from 
Eqs. (10), (111. Note that in the CO limit convergence 



takes place while the spin length j remains finite. 



III. QUANTUM CORRECTIONS TO THE 
MEAN-FIELD QPT 

The previous section showed that MF theory becomes 
exact in the CS and the CO limit, which are therefore 
identical with respect to the critical behavior of the QPT. 
The nature of the two limits is however different and they 
can be distinguished through the properties of quantum 
fluctuations around the MF ground state. 

The origin of the differences can be understood with 
a simple energy argument. In the CO limit, the energy 
scales for spin fluctuations (oc A) and oscillators fluc- 
tuations (oc 51) separate. Since A is large compared to 
the coupling constant A, which is proportional to V" 
(cf. Eq. (JsJ), spin fluctuations are suppressed in the CO 
limit. This explains partly why the QPT in the CO limit 
can occur already for finite j. In the CS limit the ra- 
tio fi/A remains constant and neither spin nor oscillator 
fluctuations are suppressed. 



A. Effective model for the CO limit 

While spin fluctuations are suppressed in the CO limit, 
oscillator fluctuations around the classical coherent state 
| a) remain energetically favorable. Their strength can 
be derived with an effective bosonic model obtained in 
perturbation theory. 

For k < 1 the MF ground state |— j')(g)|vac) is the prod- 
uct of the J z -eigenstate \—j) to the smallest eigenvalue 
— j and the bosonic vacuum |vac). For A 3> 51, the low 
energy sector of the Hilbert space consists of all states 
® |V>t>os) with a bosonic state IVw)- While the op- 
erator J z remains in the low energy sector, the operator 



J x creates a spin excitation oc \—j + 1) of energy A. 

Standard perturbation theory [14] gives the effective 
low energy model for the bosonic state IV'bos) as 



2 r 



H b 0l 



(-j\AJ z + Qa^a + ^ [(« + a ^)Jx 
-Aj + 0(a f a - -(a + a f ) 2 ) 



(14) 



Using the results from App.[A"|for the bosonic part in the 
second line, we see that the stability condition in Eq. ( A7 1 



is fulfilled only for n < 1 below the critical coupling. At 
the QPT k = 1 the number of oscillator fluctuations 



introduced through the term (a - 



it12 



diverges. 



For k > 1 we must consider spin fluctuations above 
the classical ground state \9) ® \a), which are no longer 
created by J x since 9^0. Instead, we rewrite the Hamil- 
tonian with the rotated spin operators J z — cos 9 J z — 
sm9J x , J x = sm9J z + cos9J x , and find 

H = A cos 93 z — A sin 9 (a + (v)J z + Qa^a 

+ A sin 0.4 + Acos6» (a + a r )J x . (15) 

Here the operator J x appears with the prefactor A, and 
we cannot immediately use this expression for perturba- 
tion theory for small 51/ A. 

To proceed, we shift operators J z <-> J z + j , a i— > a — a 
by their expectation values in the classical ground state, 
taken from Eqs. ([6]), ([9]), and obtain the Hamiltonian in 
the form 

H = E{9) + A(cos6» + Ksin 2 9){J Z + j) 

- Asin#(a + a f - 2a)(J z + j) 
+ £l(a f - a)(a - a) + A cos 9(a + a 1 - 2a) J x , (16) 

where the first term is the energy functional E(9) from 
Eq. Q. In this expression, a term 2a A cos 9J X has can- 
celed the problematic term Asin^J^ for the values of a, 
9 given by Eqs. @. 

Now, J x in the last line appears with a prefactor that 
is small compared to A and perturbation theory can be 
applied. Note that the spin fluctuation energy, which is 
given by the prefactor A (cos 9 + k sin 2 9) = An of J z +j, 
differs from the bare value A. The effective model for 
k > 1 is obtained as 



^bos 



+ n 



2 \k 

(at — a) (a — a) 



1 

4k 2 



(a + a 1 



2a) 



(17) 



For 



1, it coincides with 



from Eq. ( 14 ) 



The effective low energy models H£ 



describe the 

ground state of the Dicke model in the CO limit includ- 
ing oscillator fluctuations. From this model we recover 
the MF expressions for a and E(9), and thus the entire 
QPT in the CO limit. In particular the present deriva- 
tion shows that the argument given in Sec. |II C] is correct 
and not invalidated by oscillator fluctuations. 
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FIG. 5. (Color online) Spin variance Aj (left panel) and 
oscillator variance A 9 (right panel) for decreasing Q/A ap- 
proaching the CO limit, with fixed finite j — 5. The dashed 
grey curve in the right panel gives the analytical for A 9 from 
Eq. (pi). 



FIG. 6. (Color online) Entanglement entropy S in the CS 
limit (left panel) (cf. Ref. [S]) and in the CO limit (right 
panel) according to Eq. (24 1. 



IV. CRITICAL AND NON-CRITICAL 
ENTANGLEMENT 



For a quantitative analysis of the numerical data, we 
use the oscillator variance 



A 9 = (q 2 ) {qf 



(18) 



of the oscillator position q = (a + a)). Using the results 
from App. [a] the variance is obtained from H^ os as 



A 9 = 



{I-k)- 1 ' 2 

(l-l/K 2 )-^ 



if k < 1 
if k > 1 



(19) 



B. Spin and oscillator variance 

In Fig. [5] we compare the oscillator variance A q for 
small CI/ A with the spin variance Aj. We use Aj, as 
defined in App. [Cj instead of the spin variance in a fixed 
direction, for example A z = (J|) — (J z ) 2 , since it is in- 
variant under rotations. For fluctuations around the spin 
coherent state \9 = 0) that points in the z-direction it is 
identical to A z . Above the QPT, Aj accounts for the ro- 
tation of the spin axis relative to which spin fluctuations 
occur. A large Aj is an indication of significant spin 
fluctuations, while Aj — corresponds (in the present 
examples) to a spin coherent state. The vanishing of Aj 
as CI/ A — > (left panel) shows the suppression of spin 
fluctuations in the CO limit. 

For A q in the right panel we observe the growth of 
oscillator fluctuations with decreasing CI/ A. Recall that 
j = 5 is finite and small in this example, and the QPT in 
the CO limit is triggered by a macroscopic displacement 
of the classical oscillator. Oscillator fluctuations are a 
genuine quantum correction, which is independent of j 
and occurs even in the smallest non-trivial case j = 1/2. 
Approaching the CO limit, A g diverges at the QPT ac- 
cording to Eq. (19). The criticality of quantum fluctu- 
ations implies the breakdown of the classical oscillator 
limit in the vicinity of the QPT. 



In addition to the spin and oscillator variance studied 
in the previous section, corrections to the MF ground 
state arise from spin-oscillator entanglement. It can be 
measured with the entanglement entropy 



S = — Tr[plnp] 



(20) 



which is calculated with the reduced spin or oscillator 
density matrix p (both choices give the same result ac- 
cording to the Schmidt decomposition) [TSl US] . 

A simple argument would suggest a jump of S at k = 1 
from S = for the non-degenerate ground state below 
the QPT to S — In 2 for the two-fold degenerate ground 
state above the QPT. Note that we assume e = here. 

Quantum fluctuations can modify this behavior consid- 
erably, and lead to criticality of entanglement in the CS 
limit [HI [H] • In the vicinity of the QPT the entanglement 
entropy in the CS limit is given as HJ 



In II 



const. 



(21) 



such that Scs diverges at the critical coupling n — 1 with 
critical exponent 1/4. 

We show Scs in Fig. [6 (left panel). The functional form 
of Scs = Scs(k) depends on the ratio CI/ A. For CI/ A <C 
1 or CI / A 3> 1 quantum spin or oscillator fluctuations are 
energetically less favorable than for CI/ A = 1, and the 
value of S decreases away from the QPT. The criticality 
of Sqs a t the QPT and the critical exponent are however 
independent of CI/ A. 



A. Entanglement in the classical oscillator limit 

Almost trivially, the entanglement entropy S cannot 
diverge for finite j since it is bounded by S < ln(2j + 1). 
The suppression of spin fluctuations in the CO limit 
results in the much stricter condition S < In 2, inde- 
pendently of j. Large entanglement requires a sizeable 
amount of both oscillator and spin fluctuations. 

For n < 1 the ground state is a product state and 
•Scot'* < 1] = 0. This follows immediately from the fact 
used in the perturbative calculation in Sec. Ill A that the 



() 



spin part of the ground state is the single J z -eigenstate 
| — j). In contrast to the CS limit, the absence of spin 
fluctuations prevents the growth of Sco with increasing 
n. 

For k > 1 the symmetrized ground state wave function 
is given by 



\a)± 



(22) 



Note that the we assume e = here. 

The two states in the bracket are orthogonal since (ol — 
a) = in the CO limit (a diverges according to Eq. (|6j)), 
but the two spin coherent states \±9) are not orthogonal 
such that S remains strictly smaller than In 2. 

The reduced spin density matrix is 



f 



\0)(B\ + \-e)(-e\ 



(23) 



with eigenvalues fi± = 1±(9\ — 9) = l±cos 2j 9. Note that 
this expression for p s is valid also for the non-orthogonal 
states appearing here. 

The entanglement entropy obtained from p s as Sco = 
—[i- In /!_ — fi + In /i + is 



S CO [k >l]=ln2--(l- «T 2j )ln(l - K~ 2j ) 



ln(l 



(24) 



where we inserted the angle cos# = X/k according to 
Eq. ^. 

We show Sco m Fig. |6] (right panel) , where it can be 
compared to Scs- In the CO limit the entropy remains 
zero for k < 1 and increases monotonically from to In 2 
for k > 1. Note that Sqo < hi2 for all j and k. 

For finite j, it remains Sco < hi 2 even above the 
QPT because of the finite overlap (9\—9). For j —> oo, 
still strictly in the CO limit fl = 0, the overlap of the two 
spin coherent states vanishes and the curves approach a 
step function with a jump at k = 1. 



B. Quantum regime 

We show in Fig. [7] the convergence of the entanglement 
entropy S to the analytical results for the two classical 
limits. As in Figs. [3| [4] we start from the curve for j = 5, 
fi/A = 1. In the CS limit (left panel) we see how the 
divergence of S at K = 1 develops as j is increased. For 
the CO limit (right panel) we see that S remains small 
and converges to a continuous function bounded by In 2. 
Since j is finite, the limiting curve is continuous without 
the jump at k = 1 that evolves only for j — > oo. 

It should be noted that the CO limit with j > 1 and 
small entanglement differs from the regime O <C A in the 
CS limit where the entanglement entropy still diverges at 
the QPT [171 [18]. Both situations are located close to the 
origin in Fig.[l] but the first (second) situation lies closer 
to the abscissa (ordinate) than to the second axis. 





FIG. 7. (Color online) Convergence of entanglement entropy 
S towards the classical limits. Left panel: Approaching the 
CS limit with increasing j for fixed il/A = 1. Right panel: 
Approaching the CO limit with decreasing Q,/ A for fixed j = 
5. Both panel start from the curve for j = 5, Q/A — 1. The 
dashed (dot-dashed) curve gives the analytical result in the 
CS (CO) limit according to Fig. [6] and Eqs. jgj}, p4l. 



V. FAST OSCILLATOR LIMIT AND THE 
LIPKIN-MESHKOV-GLICK MODEL 

Opposed to the CO limit is the fast oscillator (FO) 
limit n/A — > oo. In this limit oscillator fluctuations 
are suppressed, while the spin fluctuations are described 
by an effective model that can be derived in perturba- 
tion theory analogously to Sec. |III A| This results in the 
Lipkin-Meshkov-Glick (LMG) model 19 known from nu- 
clear physics. For the special case j = 1/2, where the 
Dicke model reduces to the Rabi model, the FO limit can 
also be performed with a different scaling of the coupling 
constant. 



A. Derivation of the LMG model in the FO limit 

We can derive the LMG model as the effective low 
energy model for the spin part of the wave function in 
analogy to Sec. |III A] The derivation is in fact easier than 
in the CO limit, since fluctuations around the bosonic 
coherent state |a) are described by translated bosonic 
operators a) — a instead of rotated (spin) operators. The 
effective spin model, which is valid for all n, is obtained 
as the LMG model 



-Hlmg = A( J z — 



(25) 



In Fig. [8] we show the spin variance Aj (left panel) and 
the oscillator variance A q (right panel) for large fi/A. 
In reversal of the behavior in the CO limit, we see the 
suppression of oscillator fluctuations. Spin fluctuations 
remain finite and converge for f2/A — > oo to the result 
from the LMG model. 

In contrast to the CO limit no QPT and, therefore, 
no divergence of spin fluctuations occurs in the FO limit 
for finite j, simply because the initial argument against 
symmetry breaking given in the introduction applies to 
the LMG model. The QPT is recovered if additionally 
the j — > oo-limit is performed in the LMG model. Then, 
the ground state of -Hlmg is a s P m coherent state 10), 
and we recover the energy functional E{9) from Eq. (JjJ). 
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FIG. 8. (Color online) Spin variance Aj (left panel) and oscil- 
lator variance A 9 (right panel) for increasing SI/ A approach- 
ing the FO limit, with fixed large j — 50. The grey filled 
background curve in the left panel gives Aj for the LMG 
model. This figure complements Fig. [5]for the CO limit. 



Consequently, we also recover the QPT. In this sense, the 
FO and CS limit commute. 

To calculate the spin fluctuations for j — > oo, we use 
the Holstcin-Primakoff (HP) transformation [2U] 



J z = tfb-j, J x = ^J J -(tf +b)+0{r 1/2 ) (26) 

of spin operators to bosonic operators b^K For k < 1, 
when 6 = 0, the HP transformation can be applied di- 
rectly to -Hlmg and results in the bosonic model 



H<=A^b-^(b + b^-j 



(27) 



For k > 1 it is 9 ^ and spin operators must be rotated 
prior to the HP transformation, similarly to Eq. (151 in 
Sec. IIII Al We now obtain the bosonic model 



1 



^b-^{b + ^f\- J -f{n + l/n). (28) 



Comparison of these models to the effective bosonic mod 



els H^ os from Eqs. (|14|), (|17|) reveals the duality of the CO 



limit and the combined FO/CS limit, in the sense that 
the role of spin and oscillator fluctuations are reversed. 

Using the results from App. [a] for H§^ ) we find the spin 
variance Aj as (still for j — > oo) 



fl-K 



Af = 



1 

3k 2 (k 2 - 1) 



if k < 1 , 



if k > 1 



(29) 



The spin variance for large j is shown in Fig. [9] We 
observe convergence to the analytical result for j — > oo. 
Again, the QPT is accompanied by a divergence of fluc- 
tuations and a breakdown of the corresponding classical 
spin limit. 



B. FO limit for the Rabi case j - = 1/2 

The FO limit can also be treated by a unitary trans- 
formation UHW of the Dicke Hamiltonian, where U is 



FIG. 9. (Color online) Spin variance Aj for increasing j ap- 
proaching the QPT, with fixed large 17/ A = 20. The grey 
dashed curve gives the analytical result for Aj from Eq. (29 1. 



given as 



U = cxp 



- £(a f - a)J x 



A 
ft 



(30) 



The transformation U displaces oscillator states by a shift 
m x ^ that depends on the ./^-eigenvalues m x . The equiv- 
alent transformation in polaron physics is known as the 
Lang-Firsov transformation [21) . 

With the above choice for f , the interaction term (aJ + 
a)J x is eliminated through the transformation, and the 
transformed Hamiltonian reads 



UHU^ =Acosh[£(a t - a)]J z + iAsinh[£(c 
A 2 

+ Sla? a - —Jl . 
Q 



t 



l)]jy 



(31) 



So far, the transformation is an exact reformulation of 
the problem. We can now note that for $1 — > oo the pres- 
ence of the term fla^a implies that the ground state of 
UHW contains no bosonic excitations. The transformed 
ground state wave function has the form |VWn) <B> |vac). 
In the vacuum |vac) the bosonic operators from UJ Z W 
have expectation values 



(vac | cosh[£(a 1 ' — a)]|vac) = e ^ I 2 , 
(vac | sinh[£(a^ — a)]|vac) = , 

and the transformed Hamiltonian UHW reduces to 



(32) 



H 



FO 



Ae -(A 2 /n 2 )/ 2j 



(33) 



We can perform the FO limit in this model in two 
relevant ways. If we insert the coupling constant n from 
Eq. (8b as done before, we obtain the LMG model from 
Eq. (25 1. Alternatively, we can keep the parameter £ = 
X/fl of the transformation constant. Then, the prefactor 
of J 2 will diverge for Jl — > oo as we push the system into 
the strong coupling limit above the QPT. This limit is 
not interesting for general j. 

In the special Rabi case j = 1/2, however, it is 
J| = 1 /4. The first kind of FO limit is trivial for this 
model, since it results in the Hamiltonian A J z — k/4 of 
a free spin. Instead, we can perform the second FO limit 
because the divergent prefactor of J 2 now results only in 



a divergent shift of the ground state energy that can be 
dropped from the Hamiltonian. 
We then obtain the simple model 



H rC n — A J z 



(34) 



of a spin J z with renormalized frequency A = e ^ / 2 A. 
It results in the susceptibility 



Xr 



A" 



(35) 



which grows monotonically with the effective coupling 
strength £ = A/f2. The frequency renormalization is the 
sole effect of coupling to the fast oscillator. The exponen- 
tial prefactor is an example of a Franck-Condon factor, 
known from the theory of vibronic transitions or polaron 
physics. 

While the transformed model H rcn is also not very in- 
teresting by itself, with a simple ground state I — 1 / 2 ) ® 
| vac), it allows us to obtain the actual ground state of the 
Rabi model through the transformation U. It is given by 



|^) = U I-1/2) ® I vac) 



®l-0) 



(36) 



where we denote the j = ±1/2 eigenstates of J x by |— >), 
\<r-). The entanglement entropy in this state is given by 



ln2 -^(l-e- 2 « 2 )ln(l-e- 2 « 2 



i(l + e - 2 « 2 )ln(l + e - 2 « 2 ), (37) 



Sri 



(cf. the calculation for Eq. p4|). Similar to the sus- 
ceptibility, the entropy S increases monotonically with 
£ = A/Q from to In 2. 

In Fig. [lO] we contrast this behavior with the behavior 
in the CO limit (also for j — 1/2). We see that the 
entanglement entropy S is close to its maximal value In 2 
already for small x close to one. In this sense, the FO 
limit of the Rabi model is characterized by significant 
entanglement. 

Note that the relevant coupling constant scales as X/fl 
in the present FO limit but as A 2 /f2 in the CO limit. 
Previously, for the Dicke model, we had to choose the 
same coupling constant K oc A 2 /f2 for both limits, which 
leads to a duality of the FO and CO limit. For the Rabi 
model, we have defined a FO limit that is structurally 
different from the CO limit: Instead of a QPT it features 
renormalization of the spin frequency. 

The frequency renormalization is peculiar for the Rabi 
model. It is reminiscent of lattice polarons, where the 
distinction between self-trapped adiabatic polarons (cor- 
responding to the CO limit) and anti-adiabatic polarons 
(in the FO limit) involves characteristically different sig- 
natures in, e.g., the optical conductivity [22]. In fact, the 
Rabi model is equivalent to the Holstein polaron model 
restricted to two lattice sites [;23 . 
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FIG. 10. (Color online) Susceptibility \ an d entanglement 
entro py S for the Rabi case j = 1/2 in the CO limit (left pane l, 
Eqs. p4| ) and the FO limit (right panel, Eqs. {jjjj), (37)). 



VI. SUMMARY 

Two different classical limits can be defined for the 
Dicke model. The QPT in the CS spin limit j ^ 00 
has attracted much attention, and the criticality of spin- 
oscillator entanglement is understood as its characteristic 
signature. 

As pointed out here the QPT is also realized in the 
second classical limit, the CO limit SI — > 0. It should 
be noted that the QPT in the CO limit is not simply a 
special case of the CS limit: It occurs already at finite 
spin length j < 00, even at j = 1/2. 

A simple MF argument shows the equivalence of the 
QPT in the two limits if only the critical behavior is 
considered. Differences occur for quantum fluctuations 
around the MF ground state. In the CO limit, the sup- 
pression of spin fluctuations prevents significant entan- 
glement, but oscillator fluctuations are important and 
diverge at the QPT. 

The emergence of a QPT in the CO limit is a general 
feature, which can occur for any finite quantum system 
coupled to a harmonic oscillator. It does not require the 
equivalent of an j — > oo-limit. In every situation the 
QPT is accompanied by diverging oscillator fluctuations, 
while fluctuations of the finite system are suppressed. 

Also the FO limit can be performed for general sys- 
tems, although a different formulation of the limit should 
be chosen for two- level systems (the Rabi case j — 1/2) 
and systems with multiple energy levels. The duality of 
the CO and FO limit is a special feature of the Dicke 
model, where the spin can be mapped onto a bosonic 
system for j — > 00. 

The basic physical mechanisms realized in the CO and 
FO limit are typical for any finite quantum system cou- 
pled to harmonic oscillators. Oscillator fluctuations or 
entanglement with the oscillator for two-level systems in 
the FO limit can be expected to be of general importance 
in many situations. The Dicke model is an example where 
their properties can be studied in detail. One particular 
feature is the occurrence of one QPT in two classical lim- 
its, which are distinguished by entirely different quantum 
corrections to the MF ground state. 
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Appendix A: Squeezed harmonic oscillator 

The Hamiltonian of a squeezed harmonic oscillator has 
the form 



H = a f a + j3(a + a<) 2 



(Al) 



with p e R. 

The Hamiltonian can be diagonalized with a unitary 
transformation 



U = exp —(a) — a) 
for which 

U aU^ = cosh a a — sinh a a* 
With the choice 

23 



tanh2cr = 



the transformed Hamiltonian 



1 + 2/3 



H = UHU^ = y/l + 4/3a ! a + E 



(A2) 



(A3) 



(A4) 



(A5) 



acquires the form of a standard harmonic oscillator 
Hamiltonian, with ground state energy 



U 2 V H 2 1 + 4/3 



We note the stability condition 



(AT) 



For smaller /3, the original Hamiltonian is not bounded 
from below. 

Since the ground state of H is the bosonic vacuum, 
expectation values (...) in the ground state of H can be 
evaluated through transformation with U. Especially for 
the oscillator variances we find 



((a + at) 2 ) - (a + at) 2 



1 



and 



((^a) 2 ) - (<Ja) 



t„\2 = 



2/3 2 
1T4/3 



(A8) 



(A9) 



Appendix B: Some properties of coherent states 

We summarize the essential properties of oscillator and 
spin coherent states, see also Ref. |2~4") . 



Coherent states of the oscillator can be defined as the 
eigenstates of the destruction operator. 



a a) = a a) 



The coherent state can be written as 



\a) = e aa ~ a a \0) 



Expectation values are given by 

(a|a'a|a) = \a\ 2 
(a|aT + a\a) = 2Rea 



(Bl) 
(B2) 

(B3) 



2. Spin coherent states 

A (real) spin coherent state is defined as 

\0)=e i9J »\j,-j). (B4) 

It is the eigenstate of the operator cos 8 J z — sin 8 J x to 
eigenvalue —j, i.e. 



(cos8 J z - sind J x )\8) = -j\ 



(B5) 



Expectation values of spin operators in the coherent 
state are given by 



\J x \8}=j sin8, 
•\J z \8)=-jcost 



(B6) 



(A6) b y 



The overlap between two spin coherent states is given 



(8\x) = cos 



(B7) 



For j — > 00, the overlap is zero for 8 =/= \. 



Appendix C: Spin variance 

The oscillator variance A q is invariant under trans- 
lations, which modify the bosonic operators through a 
linear shift a i— > a + a. For an analogous spin variance 
we require invariance under rotations, which leads to a 
slightly more complicated definition. We restrict our- 
selves to real spin states, the generalization to arbitrary 
spin states is straightforward. 

We define the spin variance Aj as the minimal vari- 
ance of a rotated spin operator J11 = cos8J z + sin^J^, 
that is obtained through variation of the rotation angle 
8. Expansion of J| shows that Aj is the minimum of a 
quadratic form, and given by the smaller eigenvalue of 
the 2x2 matrix 



(j x )(j z ) (j 2 )-(j 2 ) 2 , 



(Cl) 
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Note that (J X J Z ) — (J Z J X ) for a real spin state. 

For a (real) spin coherent state \9), we have A,/ = 
since \ff) is obtained from rotation of the J 2 -eigenstate 
\j, —j) . Conversely, if the smaller eigenvalue A j — 0, the 
state is an eigenstate of Jn , i.e. it is a rotated eigenstate 
\j, m) of J z (the angle 9 could be deduced from the eigen- 



vectors). It need however not be a spin coherent state, 
which would require m = ±j. 

For the Dicke model in the j — > oo-limit, the spin state 
for n < 1 is invariant under the J x i— > — J x symmetry. 



The off-diagonal elements in Eq. (CI ) vanish, and Aj = 
A 2 since (Jj) is of order j. 
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